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Abstract 

Continuous-time quantum walk (CTQW) on a given graph is investigated by using 
the techniques of the spectral analysis and inverse Laplace transform of the Stieltjes 
function (Stieltjes transform of the spectral distribution) associated with the graph. It 
is shown that, the probability amplitude of observing the CTQW at a given site at time 
t is related to the inverse Laplace transformation of the Stieltjes function, namely, one 
can calculate the probability amplitudes only by taking the inverse laplace transform of 
the function iG^is), where G ll (x) is the Stieltjes function of the graph. The preference 
of this procedure is that, there is no any need to know the spectrum of the graph. 

Keywords:Continuous-time quantum walk, QD and non-QD type graphs, 
Stieltjes function, Laplace transformation 
PACs Index: 03.65.Ud 
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One of the most challenging problems in quantum computation has been the design of quantum 
algorithms which outperform their classical counterparts in meaningful tasks. Few quantum 
codes in this category have been discovered after the well-known examples by Shor and Grover 
[1, 2]. A natural way to discover new quantum algorithmic ideas is to adapt a classical one 
to the quantum model. An appealing well-studied classical idea in statistics and computer 
science is the method of random walks [3, 4]. Recently, the quantum analogue of classical 
random walks has been studied in a flurry of works [5, 6, 7, 8, 9, 10]. The works of Moore and 
Russell [9] and Kempe [10] showed faster bounds on instantaneous mixing and hitting times 
for discrete and continuous quantum walks on the hypercube (compared to the classical walk). 
The focus of this paper is on the continuous-time quantum walk that was introduced by Farhi 
and Gutmann [5]. In Refs. [11, 12, 13], the dynamical properties of CTQW on finite and 
infinite graphs have been studied by using the spectral analysis. In fact, it has be shown that 
the probability amplitudes of the CTQW on a given graph are related to the Stieltjes function 
or Stieltjes transformation of the spectral distribution associated with the adjacency matrix 
of the graph, i.e., in order to calculate the probability amplitudes one needs to evaluate the 
inverse Stieltjes transformation. The Stieltjes function has extensive applications for example 
in solid state physics and condensed matter, where it is known as Green function. It can be 
also used for calculating the two-point resistances on regular resistor networks [14] . Although, 
with any arbitrary graph one can associate a Stieltjes function (see Ref. [15]), calculation of 
the inverse Stieltjes transformation is not an easy task. In this work, we show that the 
probability amplitudes of CTQW on graphs can be evaluated by taking the inverse Laplace 
transformation of the Stieltjes function which is easier and more popular transformation with 
respect to the inverse Stieltjes transformation. In fact, the Laplace transformation is used 
extensively in various problems of pure and applied mathematics. Particularly widespread 
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and effective is its application to problems arising in the theory of operational calculus and in 
its applications, embracing the most diverse branches of science and technology. An important 
advantage of methods using the Laplace transformation lies in the possibility of compiling 
tables of direct and inverse Laplace transforms of various elementary and special functions 
commonly encountered in applications. We survey and re-derive equations for the CTQW on 
QD [12, 16] and non-QD [15] type graphs by using the Laplace transformation. 

The organization of the paper is as follows: In section 2, some preliminary facts about 
graphs and their stratifications, the quantum decomposition for the adjacency matrix of some 
particular graphs called QD graphs and Stieltjes transform of spectral measure associated 
with the graph are reviewed. In section 3, we review CTQW on an arbitrary graph and give 
a procedure for calculating the probability amplitudes of CTQWs on graphs by using the 
inverse Laplace transformation. Section 4 is devoted to some examples of CTQW on QD and 
non-QD type graphs and calculation of the corresponding probability amplitudes of the walk 
on them by employing the techniques introduced in the section 3. The paper is ended with a 
brief conclusion and an appendix containing a table of examples of CTQW on some important 
distance-regular graphs. 

2 Preliminaries 

In this section we review some preliminary facts about graphs and their stratifications, the 
quantum decomposition for the adjacency matrix of some particular graphs called QD graphs 
and Stieltjes transform of spectral measure associated with the graph. 

2.1 Graphs and their stratifications 

A graph is a pair G = (V,E), where V is a non-empty set and E is a subset of € 
j}- Elements of V and of E are called vertices and edges, respectively. Two vertices 
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i,j G V are called adjacent if G -E, and in that case we write i ~ j. A finite sequence 
i ; ii, i n & V is called a walk of length n (or of n steps) if ik-i ~ i& for all k — 1, 2, n. A 
graph is called connected if any pair of distinct vertices is connected by a walk. The degree or 
valency of a vertex x E V is defined by k(x) — \{y G V : y ~ x}\. The graph structure is fully 

represented by the adjacency matrix A defined by 

{1 if i ~ j, 
(iJeV). (2-1) 
otherwise 

Obviously, (i) A is symmetric; (ii) an element of A takes a value in {0, 1}; (iii) a diagonal 
element of A vanishes. Let h{V) denote the Hilbert space of square-summable functions on 
V, and G V} becomes a complete orthonormal basis of h{V). The adjacency matrix is 

considered as an operator acting in hiV) in such a way that 

A\t) = J2\j), ^V. (2-2) 

For i ^ j let d(i,j) be the length of the shortest walk connecting % and j. By definition 
d(i,j) — for all i G V . The graph becomes a metric space with the distance function d. 
Note that d(i,j) = 1 if and only if % ~ j. We fix a point o G V as an origin of the graph. 
Then, a natural stratification for the graph is introduced as: 

oo 

V = (J Vj(o) , V;(o) := {j G V : 9(o, j) = i}. (2-3) 

If 14(0) = happens for some k > 1, then VJ(o) = for all Z > k. With each stratum V*, we 
associate a unit vector in ^(V") defined by 

W = -L E l*>- ( 2 - 4 ) 

V ^ keVi(o) 

where, := |V^(o)| and denotes the eigenket of fc-th vertex at the stratum %. The closed 
subspace of hiV) spanned by is denoted by T(G). Since becomes a complete 

orthonormal basis of T(G), we often write 

T(G) = J2®C\<j )k }. (2-5) 
k 
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In this stratification for any connected graph G , we have 

Vi{(3) C Vi-^a) U Vi(a) U V i+1 (a), (2-6) 

for each (3 G Vi(a). Now, recall that the i-th adjacency matrix of a graph G = (V, E) is defined 

as 

f 1 if d(a,/3) =i, 
(AiUp = (a,(3EV). (2-7) 

1 otherwise 

Then, for reference state |0 O ) (|0o) = \o), with oGl^as reference vertex), we have 

M<h) = E (2-8) 

Then by using (2-4) and (2-8), we have 

AM = V^M- (2-9) 

2.2 QD type graphs 

Let A be the adjacency matrix of a graph G = (V,E). According to the stratification (2-3), 
we define three matrices A + , A_ and A as follows: for i G V k we set 

^|i> = E li>. ^-10 = E li>. A)|i>=EW- (^O) 

iey fc+ i iey fc _i jev fc 

for j ~ i. Since i 6 14 and i ~ j, then j G V^-i U14U I4+i, where we tacitly understand that 
V-i = 0. One can easily verify that 

(A + )* = A_, (A )* = A , and 

A = A + + A_ + A (2-11) 

This is called quantum decomposition of A associated with the stratification (2-3). The vector 
state corresponding to \o) = \4>o), with o G V as the fixed origin, is analogous to the vacuum 
state in Fock space. According to Ref.[16], the (A rn ) coincides with the number of m-step 
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walks starting and terminating at o, also, by lemma 2.2 of [16], if T(G) is invariant under the 
quantum components A + ,A^ and A , then there exist two sequences {u k } k X3 =1 and {a^}^ 
derived from A, such that 

A + \(j) k ) = y/uj k+l \<i) k+1 ), k>0, 

A-\<j> k ) = > 1, 

^ol^fc) = a fc+ i|0 fc ), A;>0 (2-12) 

(for more details see [11, 16]). Following Ref. [11], we will refer to the graphs with this property 
as QD type graphs and the parameters u k and a k will be called QD parameters of the graph. 
We note that in QD type graphs, the stratification is independent of the choice of reference 
state. 

2.3 Stieltjes transform of spectral measure associated with the graph 

It is well known that, for any pair (A, \4>o)) of a matrix A and a vector |0 O ), it can be assigned 
a measure /i as follows 

f i(x) = (<f> \E(x)\(l> ), (2-13) 

where E(x) = J2i \ui)(ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE(x)dx. (2-14) 
It is easy to see that, for any polynomial P(A) we have 

P(A) = J P(x)E(x)dx, (2-15) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-13) and (2-15), the expectation value of powers of adjacency matrix A over 
starting site |0 O ) can be written as 

(0o|^ m |0o) = / x m n(dx), m = 0,l,2,.... (2-16) 

JR 
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The existence of a spectral distribution satisfying (2-16) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[17], Theorem 1.2]. 

Obviously relation (2-16) implies an isomorphism from the Hilbert space of the stratification 
onto the closed linear span of the orthogonal polynomials with respect to the measure \i. From 
(4-43) and (2-9), we have for distance-regular graphs 

|0,) = ^(A)|0 O ), (2-17) 

where, P[ = -^=P% is a polynomial with real coefficients and degree i. Now, substituting (2-17) 
in (4-44), we get the following three term recursion relations for the polynomials Pj(x) 

xP' k {x) = (5 k+1 P' k+l {x) + a k P' k {x) + pkP^ix) (2-18) 

for k = 0, d - 1, with P^(x) = 1. Multiplying (2-18) by ft.../?*, we obtain 

p 1 ...p k xP k '(x) = (3 l ...(3 k+1 P' k+1 {x) + a k (3 l ...(3 k P' k {x) + fa^^x). (2-19) 

By rescaling P' k as Q k = (3i...(3 k P k , the spectral distribution /i under question is characterized 
by the property of orthonormal polynomials {Q k } defined recurrently by 

Qo(x) = l, Qi(x)=x, 

xQ k (x) = Q k+1 (x) + a k Q k (x) + PlQ^x), k > 1. (2-20) 
If such a spectral distribution is unique, the spectral distribution ji is determined by the 



identity 



GM = J*® = i = <&w = g^, (2 . 21) 

jRx-y x _ aQ ft Qd{x) f^x- Xl 



x ~ ai w 

x — a2 — 



-a 3 — ~ 

where, x\ are the roots of polynomial Qd(x). G^(z) is called the Stieltjes function (Stielt- 
jes/Hilbert transform of spectral distribution u) and polynomials {Q k ^} are defined recurrently 

as 

Q$\x) = l, QP(x)=x- ai , 
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xQ ( k \x) = Q&(aO + a k+1 Q { Z\x) + Pl^Ux), k > 1, (2-22) 
respectively. Then, from (2-21), one can see that the spectral distribution \x is written as 

d-l 

/i(x) = ^^5(x-^). (2-23) 

(=0 

The coefficients Ai appearing in (2-21) and (2-23) are calculated as 

A l= \im(x- Xl )G^x). (2-24) 
(for more details see Refs.[17, 18, 19, 20].) 

3 Continuous-time quantum walk on an arbitrary graph 

For a Given undirected graph Y with n vertices and adjacency matrix A, one can define the 
Laplacian of T as L = A — D, where D is the diagonal matrix with Djj = deg(j), the degree 
of vertex j. Classically, suppose that P(t) is a time-dependent probability distribution of a 
stochastic (particle) process on T. The classical evolution of the continuous-time walk is given 
by the Kolmogorov equation 

^ = LP®. (3-25) 

The solution to this equation, modulo some conditions, is P(t) = e tL P(0), which can be solved 
by diagonalizing the symmetric matrix L. This spectral approach requires full knowledge of 
the spectrum of L. 

A quantum analogue of the classical walk, the so-called CTQW, uses the Schrodinger 
equation in place of the Kolmogorov equation, where L is chosen as the Hamiltonian of the 
walk. This is because we can view L as the generator matrix that describes an exponential 
distribution of waiting times at each vertex. Let \<f>(t)) be a time-dependent amplitude of the 
quantum process on T. Then, the wave evolution of the quantum walk is governed by 

ihj t \<t>{t))=H\4>{t))- (3-26) 
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Assuming h — 1 for simplicity, the solution to (3-26) is given by \<f>(t)) = e~ tHt \(p(0)) which, 
again, is solvable via spectral techniques. On <i-regular graphs, we have D = |7, and since A 
and D commute, we get 

e -itH = e -it(A-\l) = jtld e -itA_ ( 3 _ 27 ) 

This introduces an irrelevant phase factor in the wave evolution. Hence we can consider 
H — A — Ai. Thus, we have 

\m = e"^|0(O)). (3-28) 

In the case of distance regular graphs, according to (4-43) the adjacency matrices are 
polynomial functions of A, hence by using (2-16) and (2-17), the matrix elements (4>i\A m | <f> ) 
for m — 0, 1, ... can be written as 

{^\A m | O > = {UP'M)A m | O ) = / x m P;(x)fi(dx). (3-29) 

JR 

By using (3-29), the probability amplitude of observing the walk at /-th stratum at time t 
can be evaluated as 

qi (t) = (^|e- M *|0 o ) = / e- iXt P;(x)u,(dx) = 1 / e^Q^^dx). (3-30) 

JR y/U)\...Ui JR 

In particular, the probability of observing the walk at starting site at time t is given by 

qo(t) = (0o|e- M *|0 o > = / e- ixt u.(dx). (3-31) 

JR 

The conservation of probability J2i=o\ (4>i\4>o{t)) | 2 = 1 follows immediately from (3-30) by 
using the completeness relation of orthogonal polynomials Pi(x). 

We notice that the formula (3-30) indicates a canonical isomorphism between the interacting 
Fock space of CTQW on distance regular graphs and the closed linear span of the orthogonal 
polynomials generated by recursion relations (4-42). This isomorphism was meant to be, a 
reformulation of CTQW (on distance regular graphs), which describes quantum states by 
polynomials (describing quantum state \<p k ) by P k (x)), and make a correspondence between 
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functions of operators (g-numbers) and functions of classical quantity (c-numbers), such as 
the correspondence between e~ tAt and e~ txt . This isomorphism is similar to the isomorphism 
between Fock space of annihilation and creation operators a, a* with space of functions of 
coherent states' parameters in quantum optics, or the isomorphism between Hilbert space of 
momentum and position operators, and spaces of function defined on phase space in Wigner 
function formalism. 

According to the result of Ref. [11], the walk has the same probability amplitudes at all 
sites belonging to the same stratum and so the evaluation of qi(t) leads to the determination 
of the amplitudes at sites belonging to the Z-th stratum Vi(o). Then, by using (2-4) we have 

qi (t) = (0*|e-^|0 o > = E = ^iHe- iM \o). (3-32) 

V Kl aeV t (o) 

Therefore, the probability amplitude of observing the walk at the site a G V/(o) at time t is 
given by 

Pa (t) = (a\e- iAt \o) = —qi(t) , for all aeVi(o). (3-33) 

3.1 Evaluation of probability amplitudes by using Laplace trans- 
form 

In this section, we give a procedure for calculating the probability amplitudes of CTQWs on 
graphs by using the Laplace transformation. To do so, first we give the definition of the Laplace 
transform. The Laplace transform of a time-dependent function /(£), denoted by f(s) = Lf(t), 
is defined as 

roc 

m = / e- st f(t)dt. (3-34) 
Jo 

The Laplace transform has the following basic properties : 
Linearity : L{af(t) + bg(t)} = af(s) + bg(s), 
Derivative : Lf'(t) = sf(s) - /(0), 



CTQW via Laplace transform 12 
Shifting: Le at f(t) =f(s-a). 

(for more information related to the Laplace transform see [32], [22] and [23]). 

By taking the Laplace transform of the amplitudes qi(t) in (3-30) and using the equation 
(2-23), one can obtain 

I roc r I roo 

qi{s) = - / e~ st / e'^Q^^dt = ]T A^x,) / e^'^dt = 

t/ui—Ui Jo Jr xM.-.c^ ^ Jo 

1 ^ AjQi(xi) = i ^ A iQi( x i) ( 3 35) 



L - J S + iXi y/0Ji...UJi Y is - x i 

Therefore, although it is not an easy work in the most cases, the probability amplitudes qi(t) 
of the walk can be obtained by taking the inverse Laplace transform of (3-35), but we have no 
need to do so. Instead, we calculate only qo(s) as follows 

roo r -A- 

q (s) = / e - st / e- ixt n{dx)dt = iV — = iGJis). (3-36) 

Jo Jr ~ is — Xi 

and obtain the probability amplitude of observing the walk at starting site at time t as 

q (t)=iL- 1 (G li (is)). (3-37) 

From the formulas (3-30) and (3-31), it can be easily seen that 

ft(t ) = W - )fc(t ) = _^_g l(i (3-38, 

In fact, we need only to know the probability amplitude of observing the walk at starting site at 
time t and the polynomials Qi(x) which are obtained via the recursion relations (2-20). Then, 
by using the equations (3-37) and (3-38), the probability amplitudes qi(t) can be evaluated. 

It could be noticed that, in this approach we do not any need to take the inverse Stieltjes 
transform to obtain the spectral measure //; Instead we take the inverse Laplace transform of 
the Stieltjes function G^{is) which is more popular and convenient (there are many handbooks 
and tables of Laplase transforms [32], [22], [23] and computer programmings for this purpose) 
in comparison with the inverse Stieltjes transform. In the following section, we calculate the 
probability amplitudes for CTQW on some important QD graphs, i.e., distance-regular graphs 
and some non-QD type graphs, explicitly. 
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4 Examples 

In this section we consider CTQW on some examples of graphs and calculate the corresponding 
probability amplitudes of the walk on them by employing the techniques introduced in the 
previous section. 

4.1 Examples of distance-regular graphs 

In this subsection first we consider some set of important QD graphs called distance-regular 
graphs. To do so, we recall the definitions and properties related to distance-regular graphs: 
Definition. An undirected connected graph G = (V, E) is called distance-regular graph 
(DRG) with diameter d if it satisfies the following distance-regularity condition: 
For all h, i,j G {0, 1, d}, and a, (3 with d(a, f3) = h, the number 

V% =| {7 e V : d(a, (3) = i and 9( 7 , (3) = j} | (4-39) 

is constant in that it depends only on h, i,j but does not depend on the choice of a and (3. 
Then, := \Vi(a) n Vj{(3)\ for all a, (3 E V with d(a,P) = h. This number is called the 
intersection number. In a distance regular graph, p % - x = (for i ^ 0, j is not {i — + 1}). 
The intersection numbers of the graph are defined as 

Oi=Pii, b i =p\ +lil , Ci=p\_ 1A . (4-40) 
These intersection numbers and the valencies Ki satisfy the following obvious conditions 

Q>i ~\~ h% + Cj = K, Ki—ibi—i = KjCii i — 1, d, 

k = ci = 1, b = k x = k, (co = b d = 0). (4-41) 

Thus, all parameters of the graph can be obtained from the intersection array {b , 6^-1! c ii ■■■■> c d}- 
It could be noticed that, for adjacency matrices of a distance regular graph, we have 

A x Ai = + ciiAi + c i+1 A i+1 , for i = 1, 2, d - 1, 
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A x A d = b d -\A d -\ + a d A d . 



(4-42) 



Using the recursion relations (4-42), one can show that Ai is a polynomial in A x of degree i, 
i.e., we have 

A i = P i (A 1 ), i = l,2,...,d, (4-43) 

and conversely A\ can be written as a linear combination of I, Ai, A d (for more details see 
for example [11]). 

It should be noticed that, for distance-regular graphs, the unit vectors \<pi) for i — 0, 1, d 
defined as in (2-4), satisfy the following three-term recursion relations 



A\(j>i) = p i+1 \(/> i+1 ) + ailfc) + 



(4-44) 



where, the coefficients a, and (3i are defined as 



a k = k - b k - Cfc, u k = 01 = b k -ic k , k = l,...,d, 



(4-45) 



i.e., in the basis of unit vectors = 0, l,...,d}, the adjacency matrix A is projected to 

the following symmetric tridiagonal form: 

/ _ a n n \ 



A 



a (3i 
Pi «i Pi 

Pi «3 Pz 















(4-46) 



. p d -i a d -i p d 
V U ... p d a d ) 

In the following, we consider CTQW on some examples of distance-regular graphs in details, 
where in the Appendix A, we give a table of important distance-regular graphs together with 
the corresponding probability amplitude qo(t). 
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4.1.1 Complete graph K n 

The complete graph K n is the simplest example of distance-regular graphs. This graph has n 
vertices with n{n— 1)/2 edges, the degree of each vertex is k = n—1 also the graph has diameter 
d—1. The intersection array of the graph is {60; Ci} = {n — 1; 1}. Then, the graph has only- 
two strata V (a) = a and Vi(a) = {(5 : (3 7^ a} with QD parameters {o^; cui} = {n — 2, n — 1}. 
By using (2-21), the Stieltjes function is calculated as 

From (3-36) and (4-47), we obtain 

- s + ^( n — 2) 1 1 n — 1 

* (S) = 8 2 + i 8 (n-2)+n-l = A + ^-l) + ~ } ' (4 " 48) 

Therefore, by taking the inverse Laplace transform of (4-78), we obtain the probability ampli- 
tude of observing the walk at starting site at time t as follows 

q (t) = L-^qois)) = I( e - < (»- 1 )* + (n - l)e«). (4-49) 
n 

For calculating the probability amplitude qi(t), first recall that k± — n — 1 and = ^7=bj> 
then by using (3-38) we obtain 

qi (t) = iL(f) = ^EI (e -i(n-l)* _ ^ (4 . 50) 

Vn — 1 at n 

4.1.2 Strongly regular graphs 

For strongly regular graphs we have three strata and two kinds of two-point resistances R a/3 (i) 
and R a/3 (2). The QD parameters of the graph with parameters (v, k, A, /i) are given by 

«i = A, a 2 = k — /i] lo\ — k, 0J2 = — A — 1). (4-51) 

Then, by using (2-21) and (3-36), one can obtain 

~ ( \ = s 2 + i(K + \ - fi)s - k{\ - n) - fi 

qo{S) s 3 + i(k + A - fi)s 2 + (k(h - A + 1) - fj)s + ik(k -h)' { } 
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Therefore, the probability amplitude qo(t) can be obtained as follows 

q (t) = A ie - ZKt + A 2 e-2( A -' 1 +V( A -^) 2 - 4 ^- "))* + A3e -^-»W( x -») 2 -^-K))\ (4-53) 

where, 



k 2 — k(A — //) + (// — «)' 



-KyJ(X - /l) 2 - 4(/i - /c) + k(A - /i) + 2/c 
= 



(A - - 2k) v /(A - if) 2 - 4(/j - k) + (A - if) 2 - 4(/j - k) 

kJ(X - it) 2 - 4(it - k) + k(A - //) + 2k 

A 3 = v . — . (4-54) 

(-A + n + 2/t) v /(A - if) 2 - 4(u, - k) + (A - if) 2 - 4(u, - «) 

In the remaining part of this example, we study CTQW on the following two well-known 
strongly regular graphs: 

A. Petersen graph 

The Petersen graph [27] is a strongly regular graph with parameters (v, k, A, //) = (10, 3, 0, 1), 
the intersection array {6o, b\\ c±, c 2 } = {3, 2; 1, 1} and QD parameters {a±, a 2 ; u>i, ou 2 } = {0, 2; 3, 2}. 
Then, by using (4-54), we obtain 

A 1 = —, A 2 = -, A 3 = -. (4-55) 
10 2' 3 5 V ; 

By substituting (4-55) in (4-53), we obtain the probability amplitude of observing the walk at 
starting site at time t as follows 

qo (t) = ^(5e-** + Ae m + e~ m ). (4-56) 

Then the probability amplitudes qi(t) and q 2 (t) are easily calculated as 

*<*> = 75'-^ - 3) * (i) = 76 ( - e "" + r* + r 3,,) - (4 - 57) 

B. Normal subgroup scheme 

Definition The partition P = {P , Pi,..., Pd} of a finite group G is called a blueprint if 
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0) = {e} 

(ii) for i=l,2,...,d if g G Pi then g~ x G P { 

(iii) the set of relations P; = {(a, j3) G G^G^ 1 ^ G p} on G form an association scheme [2 7]. 
The set of real conjugacy classes given in Appendix A of Ref. [11] is an example of blueprint 
on G. Also, one can show that in the regular representation, the class sums Pj for i — 0, 1, d 
defined as 

P i =Y,lzCG, i = 0,l,...,d, (4-58) 

are the adjacency matrices of a blueprint scheme. 

In Ref. [11], it has been shown that, if H be a normal subgroup of G, the following blueprint 
classes 

P = {e}, P 1 = G-{H}, P 2 = Pf-{e}, (4-59) 

define a strongly regular graph with parameters (v, k, A, /i) — (g, g — h, g — 2h, g — h) and the 
following intersection array 

{b , bi, a, c 2 } = {g - h, h - 1; 1, g - h}, (4-60) 

where, g := \G\ and h := \H\. It is interesting to note that in normal subgroup scheme, 
the intersections numbers and other parameters depend only on the cardinalities of the group 
and its normal subgroup. By using (4-45), the QD parameters are given by {o?i, oli\ uj±, co^} = 
{g-2h,0;g-h, (g-h)(h-l)}. 

As an example, we consider the dihedral group G = D 2m , where its normal subgroup is 
H = Z m . Therefore, the blueprint classes are given by 

P = {e}, P 1 ^{b,ab,a 2 b,...,a m - 1 b}, P 2 = {a,a 2 , ...,a^}, (4-61) 

which form a strongly regular graph with parameters (2m, m, 0, m) and the following intersec- 
tion numbers and QD parameters 

{b ,b 1 ;c 1 ,c 2 } — {m, m - 1; l,m}; {a 1 , a 2 ; u 2 } — {0, 0; m, m(m - 1)}. (4-62) 
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Then, by using (4-54), we obtain 

Now, by substituting (4-63) in (4-53), we obtain the probability amplitude of observing the 
walk at starting site at time t as follows 

q (t) = — (m - 1 + cos mi). (4-64) 

The other probability amplitudes are calculated as 

1 d 2 

q 2 (t) = —=( — — - my/m)q (t) = (y/m - 1) cosmt — (m - 1). (4-65) 

\J TYl ut 

4.1.3 Cycle graph C 2m 

The cycle graph or cycle with n vertices is denoted by C n with k — 2. We consider n = 2m 
(the case n = 2m + 1 can be considered similarly). The intersection array is given by 

{b , 6 m _i; ci, c m } = {2, 1, 1, 1; 1, .., 1, 2} (4-66) 

Then, by using (4-45), the QD parameters are given by 

cui = 0, i = 0, 1, m; ui = u m = 2, — i — 2,...,m—l, (4-67) 

Therefore, by using (2-21) and (3-36), one can obtain 

• rpf(is\ 

* w - *m- « 4 - 68) 

Then, the probability amplitude q (t) is obtained as 

i 2n-l 

q (t) = -(cos2t + e- 2ltcos ^). (4-69) 

n 1=1, tyn 
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From (3-38), one can calculate 

= _L| (90(i) ) = _^ (sill2i + ! J£ cos * £ --^). (4 . 7 0) 

The amplitudes qi(t), for / > 1 can be calculated similarly, where the results thus obtained are 
in agreement with those of Ref.[25]. 

It could be noticed that, in the limit of the large n, the cycle graph tend to the infinite line 
graph Z and the Stieltjes function reads as 

G,(x) = -i4=. (4-71) 
yx l — 4 

Therefore, by using (3-36) we obtain q (t) as follows 

qo (t) = L-\iG,{is)) = J (2t), (4-72) 
where, the J is Bessel function. From (3-38), we can calculate 

qi{t) = i/Ut {U2t)) = 

= 7f(-^ - 2)(^o(2t)) = - V2J 2 (2t). (4-73) 
By using (4-87), one can deduce that 

qi (t) = V^i-ip^t), (4-74) 

where the results are in agreement with those of Ref.[26]. 

4.1.4 Johnson graphs 

Let n > 2 and d < n/2. The Johnson graph J(n,d) has all <i-element subsets of {1,2, ...,n} 
such that two rf-element subsets are adjacent if their intersection has size d — 1. Two ci-element 
subsets are then at distance % if and only if they have exactly d — i elements in common. The 
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Johnson graph J(n, d) has v = d ^[ d ^ vertices, diameter d and the valency k = d(n — d). Its 
intersection array is given by 

bi = (d-i)(n-d-i), 0<i<d-l; Ci = i 2 , \<i<d. (4-75) 

Then, by using (4-45) the QD parameters are given by 

a k = k{n-2k), uo k = k 2 {d - k + l)(n - d - k + 1). (4-76) 

For d = 2, the Stieltjes function is calculated as 

G " (a ° = x 2 -(n-2)x-2(n-2)- ( } 

Then by using (3-36) and (4-47), we obtain 



. / x = s + t(ra-2) = 1 1 ~ y^i ! 1 + VT^I 

9 °' S ' ~~ S 2 + is(n - 2) - 2(n - 2) ~~ 2^ , ., n-2+y/(n-2)(n+6) ^ , v n-2- v /(n-2)(n+6) /' 

S + H 2 J S + H 2 / 

(4-78) 

Then, the probability amplitude qo(t) can be obtained as follows 



. M , y (n-2)(n + 6) , ^(n-2)(n + 6) 
gb(t) = L ^G^is)) = e 2 *{cos(^ )f + l y__ sm (V _ 

(4-79) 



4.2 Examples of non-QD type graphs 
4.2.1 Tchebichef graphs 

By choosing Tchebichef polynomials of the first kind (the second kind )with scaling factor 
l/2 m as orthogonal polynomials appearing in the recursion relations (2-20), i.e., Q n (x) = 
2^- 1 ) n + 1 T n (x/2 m ) ( y 2 ( - m ~ 1 ^ n U n ( y x/2 m )), one can obtain a class of finite and infinite graphs of 
Tchebichef type, with QD parameters ui = 2 2 ( m " 1 ) +1 , u k = 2 2( - m ~ l \ k = 2,3, and a k = 0, 
for k = 1,2, 3, ... {a k = 0;u k = 2 2{ - m ~ l \ k = 1,2, ...). Then, we obtain 

* w = SiuM '^mJi'' (4 " 80) 
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Therefore, the probability amplitude qo (t) can be obtained as follows 

go(t) = -Ee- l2micosi ^ (4-81) 



1=0 

for the first kind and 



for the second kind. 



4.2.2 Finite path graph P n 

For m = 1 in (4-80) and Tchebishef polynomials of the second kind we obtain finite path graph 
P n = {0, 1,2, ...,n — 1}, where it is a n-vertex graph with n — 1 edges all on a single open 
path [11]. For this graph, the stratification depends on the choice of the starting site of walk. 
If we choose the first vertex as the starting site of the walk, the QD parameters are given by 
uji — 1, ctj = 0, i — 1, ...,n — 1. Then by using (4-82) we have 

q Q (t) = Vsin 2 (-^-)e^ cos ^). (4-83) 

From (3-38) and the recursion relations P( +1 (x) = xP[(x) — P/_ 1 (x), one can obtain the prob- 
ability amplitude of observing the walk at Z-th stratum at time t as follows 

where the results thus obtained are in agreement with those of Ref. [26]. 

In the limit of the large n, one can calculate the Stieltjes function as follows 

X 

then we obtain the following closed form for the Stieltjes function 



x _ J x 2 - 4 

G,{x) = \ . (4-86) 



CTQW via Laplace transform 22 
Therefore, by using (3-36) we have 

«.) = ^H^i (4-87) 
and then the probability amplitude of observing the walk at starting site at time t is given by 

qo(t) = J (2t) + J 2 (2t). (4-88) 
Now, by using (3-38), the probability amplitudes qi(t) are obtained as follows 

q l (t)=i\j l (2t) + J l+2 (2t)). (4-89) 

4.2.3 The graphs G n 

For m — | in (4-80) and Tchebishef polynomials of the second kind we obtain a sequence 
of graphs denoted by G n (for more details see Ref.[30]). The number of vertices in G n is 
2™ +1 + 2™ — 2. In general, G n consists of two balanced binary trees of depth n with the 2™, n-th 
level vertices of the two trees pairwise identified. For the quantum walk on G n , we assume 
that the starting site of the walk is the root of a tree and calculate the probability amplitude 
of the presence of the walk at the other vertices as a function of time. Clearly, the graph G n 
has (2n + 1) strata, where the j-th stratum consists of 2? vertices for j = 1,2, ...,n + 1 and 
22n+i-i f or j _ n i ^ 2 n + 1. The QD parameters are u)i — 2, at = 0, % — 1, In. 
By using (4-82) we obtain 

*W~|^(^)^- ( -'. (4-90) 

From (3-38) and the recursion relations P[ +X {x) = xP((x) — 2P[_ 1 (x), one can obtain the 
probability amplitude of observing the walk at /-th stratum at time t as follows 

^~4r ( ^ )sH ^ )e ~' 2 ^ & ' 1> - L (4 " 91) 

In the limit of the large n, one can calculate the Stieltjes function as follows 
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then we obtain the following closed form for the Stieltjes function 

G,(x) = x ~^f^ . (4-93) 



Therefore, by using (3-36) we have 



and then the probability amplitude of observing the walk at starting site at time t is given by 

q (t) = J (2V2t) + J 2 (2V2t). (4-95) 
Now, by using (3-38), the probability amplitudes qi(t) are obtained as follows 

qi (t) = i l (J t (2V2t) + J l+2 (2V2t)). (4-96) 

4.2.4 Finite path graph with the second vertex as starting site of the walk 

Now, we study an example of non-QD type graphs such that as the authors have shown in 
Ref. [15], one can give the three term recursion to the graph by using the Krylov-subspace 
Lanczos algorithm applied to the chosen reference state and the adjacency matrix of the graph 
and calculate the Stieltjes function (for more details see [15]). For instance, we consider the 
finite path graph and choose the second vertex of the graph as the starting site of the walk, 
then the graph does not satisfy three term recursion relations, i.e., the adjacency matrix has 
not tridiagonal form. 

In Ref. [15], it has be shown that the QD parameters for P n are given by q.{ = for 
i — 0, 1, 2k — 1 and 

w 2i = ^4, W 2 i-l = ' t j^, U)2k-1 — \\ l = l,...,k — l. 

for even values of n, where for odd values of n we have 
^2i = -pr, i = l,...,fc-l. 
wm-i = lJ r, i = l,...,k, 
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Substituting these QD parameters in (2-20) and (2-22) and using (2-21), the Stieltjes function 
is obtained as 

xU n .. 2 (x/2) 

MJ U n (x/2) (49 ° 
where, C/ n 's are Tchebichef polynomials of the second kind. Therefore, by using (3-36) we have 



A -sU n - 2 (is/2) 

q ° {s) = U n (is/2) ■ (4 " 98) 
Then, the probability amplitude of the walk at starting site (the second vertex) at time t is 
given by 

1 n 9/ -7T 

Tl ~\~ 1 Ti 1 

again one can calculate the other probability amplitudes by using the Eq.(3-38). 



5 Conclusion 

CTQW on graphs was investigated by using the techniques based on spectral analysis of the 
graph and inverse Laplace transformation of the Stieltjes function associated with the graph. 
It was shown that the probability amplitudes of observing the walk at a given site at time t can 
be evaluated only by calculating the inverse Laplace transformation of the function iG^{is) 
(without any knowledge about the spectrum of the graph), where G^{x) is the Stieltjes func- 
tion associated with the graph. 



Appendix A 

In this appendix, we give the Laplace transform of the probability amplitude of observing the 
walk at starting site at time t, qo(s), for some important distance- regular graphs with n < 70. 



CTQW via Laplace transform 



25 



The graph with Ref. 


Intersection array 


9o( s ) 


90 (*) 


Icosahcdron[28] 


{5, 2, 1; 1, 2, 5} 


-i(is 3 -4s 2 +5is-10) 
s 4 +4is 3 + 10s 2 +20is + 25 


Y^(5e U +e^ 5it +6cosx/5t) 


L(Petcrscn)[28] 


{4, 2, 1; 1, 1, 4} 


-i(is 3 -3s 2 +4is-4) 
s 4 +3is 3 +8s 2 + 12is + 16 


Tj^e" + e" 4it + 10cos2t) 


Pappus, 3-covcri^3 3 [28] 


{3, 2, 2, 1; 1, 1, 2, 3} 


s 4 + 9s 2 +6 
s(s 4 + 12s 2 +27) 


^ (cos 3t + cos \/3t + 2) 


IG(AG(2, 4) \ pc)[30] 


{4, 3, 3, 1; 1, 1, 3, 4} 


s 4 + 16s 2 + 12 

sfs4ion q 2if,4^ 
s^ T^" 1 T U4 i 


i (cos 4i + 12 cos 2t + 3) 


3- covering [29] 


{9, 8, 6, 1; 1, 3, 8, 9} 


s 4 + 81s 2 +216 
sfs 4 + 90s 2 +7291 


J=r (cos 9i + 18 cos 3t + 8) 


Odd(4)[32] 


{4, 2, 1; 1, 1, 4} 


i(s 4 +3is 3 + 13s 2 +15is + 20) 
is 5 -3s 4 + 18is 3 -30s 2 +65is- 75 


Y5(4e iJ + e- 4it + 10cos2t) 


SflG\sprcad[34] 


{9, 6, 1; 1, 2, 9} 


-i(is 3 -8s 2 +9is-18) 
s 4 +8is 3 + 18s 2 + 72is + 8 


3^(96** + e - 9it + 30cos3t) 


3-covcrK 6i6 [29] 


{6, 5,4, 1; 1,2, 5, 6} 


-i(is 3 -16s 2 + 17is-136) 
s 4 4-16zs 3 +3 4 s 2 +272is + 289 


^ (2 cos 6t + 24 cos Vfit + 10) 


Hadamard graph [29] 


{12, 11, 6, 1; 1, 6, 11, 12} 


s 4 + 144s 2 +79 2 
s(s 4 + 156s 2 + 1728) 


Jj(cos 12t + 12 cos 2VEt + 8) 


IG(AG(2,5) \pc)[30] 


{5, 4, 4, 1; 1, 1, 4, 5} 


s 4 + 25s 2 +20 
s(s 4 + 30s 2 + 125) 


^(cosSi + 20 cos \/3i + 11) 


Hadamard graph[31] 


{8,7,4,1; 1,4, 7, 8} 


s 4 + 64s 2 +224 
s(s 4 + 72s 2 +51 2 ) 


Jj(2cos8t + 16cos2\/2t + 14) 


Dcsargucs[28] 


{3, 2, 2, 1, 1; 1, 1, 2, 2, 3} 


s(s 4 + lls 2 +22) 
s u + 14s 4 + 49s 2 4-36 
-i(is 3 -6s 2 +7is-14) 
s 4 +6is 3 + 14s 2 +42is + 49 
-i(is 3 -6s 2 +3is-24) 
s(s 3 + 6is 2 +9s + 54i) 
i(s 4 +2is 3 +5s 2 + eis + 2 ) 
is 5 -2s 4 +8is 3 — 12s 2 + llis-6 
-i(is 3 -12s 2 + 13is-5 2 ) 
s 4 + 12/- :i 1 2li.'- 2 I ir,(j/w + 169 


^ (cos 3t + 4 cos 2t + 10 cos t) 


Klein[28] 


{7, 4,1; 1,2, 7} 


Jj(7e" + e - 7it + 16cos\/7i) 


_ff(3,3)[28] 


{6, 4, 2; 1,2, 3} 


J-( e - 6 " + 8e 3 " + 6e- 3 " + 12) 


coxeter[28] 


{3, 2, 2, 1; 1, 1, 1, 2} 


5g(19e" + 8e~ 2 " + 12 cos s/2t) 


Mathon(Cycl(13, 3))[35] 


{13, 8, 1; 1,4, 13} 


i(13e" + e - 13 « +28COSV/T3*) 


Taylor(P(17))[29] 


{17, 8, 1; 1, 8, 17} 


-i(is 3 -3s 2 +4is-4) 
s 4 +3is 3 ■ 8.^ 2 | 12is + 16 


Jg(17e" + e - 17it + 18 cos vT7t) 


Taylor(SRG(25, 12))[29] 


{25, 12, 1; 1, 12, 25} 


-i(is 3 -24s 2 + 25is -300) 
s 4 + 2 1j : no.-- 2 (i I H t ; .-- + 625 


^(25e" + e - 25 « +26cos5t) 


lvlathon(Oycl(16, 3JJ[3oJ 


{16, 10, 1; 1, 5, 16} 


-i(is 3 - 15s 2 + 16is -80) 
s 4 4-15is 3 +32s 2 +240is + 256 

-i(is 3 -10s 2 + llis-22) 
s 4 + 10is 3 +22s 2 + 110is + 121 


(16e" + e~ leit + 34 cos it) 


Mathon(Cycl(ll, 5))[35] 


{11, 8, 1; 1, 2, 11} 


jL(lle" Xe- 11 " + 48 008^*) 


Mathon(Cycl(19, 3))[35] 


{19, 12, 1; 1, 6, 19} 


-i(is 3 -18s 2 + 19is-114) 
s 4 4-18is 3 +38s 2 +3 42 is + 361 


^(19e" +e- 19it +40COSV/T9*) 


Taylor(SflG(29, 14))[29] 


{29, 14, 1; 1, 14, 29} 


-i(is 3 -28s 2 +29is-406) 
s 4 4-28is 3 +58s 2 +812is + 8 4 l 


gL(29e" + e - 29it +30cosv^9t) 


Taylor(P(13))[29] 


{13, 6, 1; 1, 6, 13} 


-i(is 3 -12s 2 + 13is-78) 
s 4 4-12is 3 +26s 2 + 156is + 169 


^(13e" + e - 13 " + 14cos vT3t) 


GQ(2,4)\sprcad[28] 


{8, 6, 1; 1, 3, 8} 


-i(is 3 -5s 2 +22is-8) 
s 4 +5is 3 +30s 2 +40is + 64 


i(8e" + e - 8it + 12e- 2 " + 6e 4it ) 


Doro 


{12, 10, 3; 1, 3, 8} 


-i(is 3 -1 Is 2 +20is- 120) 


i( e - 12it + 17e- 4 " + 16e 6 " + 34) 




s(s 3 + llis 2 +32s + 240i) 


Locally Petersen 


{10, 6, 4; 1, 2, 5} 


-i(is 3 -12s 2 -15is-60) 

S{S -\- LZZS — OS-f- lOU! ) 


J-( e - ltHt + 13e- 5it + 25e 3it + 26) 


Taylor(GQ(2, 2)) [29] 


{15, 8, 1; 1, 8, 15} 


is 3 -12s62+43is-90 
s 4 + 12!s^+58s 2 + 180is + 225 
-i(is 3 -16s 2 -13is-120) 
s 4 + 16zs 3 +2s 2 '+240is + 225 


Jj(15e" + 6e 5it + lOe" 3 ** + e - lsit ) 


Taylor(T(6))[29] 


{15, 6, 1; 1, 6, 15} 


^(15e" + 10e 3it + 6e- 5it + e - lsit ) 


Gosset,Tayl(Schlfifli)[29] 


{27, 10, 1; 1, 10, 27} 


-i(is 3 -32s 2 -129is-432) 

= 4 1 OO ; . .i 1 I)-) .- IJaXTOO 

s Tjiis — j_ u z s -j-on^is-]- f 


i(27e it + e" 27 " + 7e~ 9 " + 21e 3 ") 


Taylor(Co-Schlafli)[29] 


{27, 16, 1; 1, 16, 27} 


-i(is 3 -20s 2 + 183is-270) 
s 4 20; | 210.- 2 | r,.H)i .s + 729 


i(27e" + e- 27 " + 7e 9it + 21 e - 3it ) 


Gfl"(2, 2)[29] 


{6, 4, 4; 1,1, 3} 


-i(is 3 -5s 2 +9is-21) 
s 4 +5i.^ :i 1- ir,s 2 |--15is + 54 


Jj^e" + e- eit + 14e 3it + 21 e - 3it ) 


ff(3,4),Doob[29] 


{9, 6, 3; 1, 2, 3} 


-i(is 3 -12s-23is-42) 
s 4 + 12is 3 — 14s 2 + 132is — 135 


5^(276"" + 27e 3 " + 9e~ 5 " + e" 9it ) 


Wells[30] 


{5,4, 1,1; 1,1, 4, 5} 


i(s 4 +3is 3 + 13s 2 +15is + 20) 
is 5 — 3s 4 + 18is 3 — 30s 2 +65is — 75 


^ (10e _it + e~ 5it + 5e 3it + 16 cos V5t) 


Gfl"(2, 1)[28] 


{4, 2, 2; 1, 1, 2} 


-i(is 3 -4s 2 +is-6) 
s 4 + 4t.v' ■ r>-- 2 | I8is-f8 
-z(is 3 -8s 2 -llis-8) 
s 4 +8is 3 - 5s 2 +44is-12 


Jj(e- 4it + 8e 2it + 12( e -" + cos V2t)) 


Gfl"(3, 1)[36] 


{6, 3, 3; 1, 1, 2} 


■Jj( e - 6it + 27e 2 " + 24e- 2it + 24 cos VZt) 


Do decahedron [28] 


{3, 2, 1, 1, 1; 1, 1, 1, 2, 3} 


-z(is 5 -2s 4 + 7is 3 -19s 2 + 10is-6) 
s(s 5 +2is 4 + 10s 3 + 16is 2 +25s + 30i) 


± (5e _it + 4e 2 " + e" 3 " + 6 cos VEt + 4) 


Perkel[37] 


{6, 5, 2; 1, 1, 3} 


-i(is 3 -6s 2 +2is-15) 
s 4 + 8s 2 +6is 3 +51is-18 


^( e - 6it + 20e 3 " + 36e- 3 "/ 2 + 36 cos V5/2i) 


GO(2, 1)[30] 


{4, 2, 2, 2; 1, 1, 1, 2} 


i(s 4 +5is 3 -s 2 + 13is-2) 
is 5 -5s 4 + 3is 3 -29s 2 +2is-24 


■±(9e u + 10e-" + 16e 2,t + 9 e - 3it + e" 4 ") 


3-coverGQ(2, 2) [30] 


{6, 4, 2, 1; 1, 1, 4, 6} 


i(s 4 +5is 3 + lls 2 + 33is + 6) 
is 5 -5s 4 + 17is 3 -57s 2 +72is-108 


±(9e- it + 18e 2it + 5e 3it + 12 e - 3 " + e" 6 ") 


J(8,4)[29] 


{16, 9, 4, 1; 1, 4, 9, 16} 


i(s 4 +20is 3 -44s 2 +368is-192) 
is 5 -20^ i 2Ho :! r,<>2.-- 12S/.S-2048 


^L_ (e -16it + 7e -8»t + 2S e 2it + 20e" 2it + 14e 4 " ) 
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